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■ Abstract. The free energy of a closed 3- manifold is a 2-parameter formal power series 

which encodes the pcrturbative Chern-Simons invariant (also known as the LMO invariant) 
of a closed 3- manifold with gauge group U (N) for arbitrary N. We prove that the free energy 
of an arbitrary closed 3-manifold is uniformly Gevrey-1. As a corollary, it follows that the 
genus g part of the free energy is convergent in a neighborhood of zero, independent of the 
genus. Our results follow from an estimate of the LMO invariant, in a particular gauge, and 
from recent results of Bender-Gao-Richmond on the asymptotics of the number of rooted 
maps for arbitrary genus. We illustrate our results with an explicit formula for the free 



energy of a Lens space. In addition, using the Painlcve differential equation, we obtain an 

-d 



(-h ! asymptotic expansion for the number of cubic graphs to all orders, stengthening the results 

of Bender-Gao-Richmond. 
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G\ ■ 1 Introduction 

O 

^ ■ 1.1 The free energy of a closed 3-manifold 

The free energy Fm(t, H) of a closed 3-manifold M (defined in equation (|1.7p ) is a 2-parameter 
^ \ formal power series with rational coefficients in two variables r and h 2 

F m (t, h) = Y^ ^(r)^- 2 G ^ 2 Q[[r, ft 2 ]]. (1.1) 

The variable h plays the role of Planck's constant, and the variable r = Nh is the product 
of Planck's constant with N, the size of the gauge group U(N). The free energy encodes the 
perturbative Chern-Simons invariant of M along the trivial flat connection (also known as the 
LMO invariant of [U]) with gauge group U(N) for arbitrary N. Fm^t) (resp. Fm,o{t)) is often 
called the genus ^-contribution (resp. planar limit) to the free energy. Perturbative Quantum 
Field Theory for a gauge group of fixed size N typically leads to factorially divergent formal 
power series partly because there are factorially many Feynman diagrams, as was explained for 
example in [20]. It is a fundamental question how to give rigorous analytic meaning to these 
factorially divergent series and how to numerically evaluate them for example for QCD, and 
compare with well-tested experimental data. With this in mind, 't Hooft observed that when 
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the size N of the gauge is arbitrary, the Feynman diagrams of perturbative gauge theory organize 
themselves in ribbon graphs, i.e., abstract connected oriented surfaces of some genus g with some 
nonzero number of boundary components; see [35]. Correspondingly, the free energy Fm{t,H) 
becomes a sum of power series FM,g(j) of a single variable that carries the contribution of the 
connected graphs of genus g. 't Hooft suggested that 

Conjecture 1.1. There exists a disk Dm that contains zero such that for all g, Fm >9 (t) is 
a power series analytic in Dm- 

If the above conjecture holds, then Fm, 9 (t) can be evaluated near r = by a convergent series. 
The above conjecture can be formulated for any gauge theory, 't Hooft's conjecture has been 
verified in some simple quantum-field theoretical models with U(N) symmetry, such as matrix 
models (QFT in zero dimensions) and matrix quantum mechanics (QFT in one dimension); 
see p2J [101 ESI H3 [55]. Evidence for its validity in N = 4 super Yang-Mills theory has also 
appeared recently; see [9[ [8]. 

1.2 A free energy of a knot 

In the world of Quantum Topology, the free energy of a 3-manifold is structurally similar to two 
other invariants of knots: namely the HOMFLY polynomial Hx(t, q) of a knot K and the Colored 
Jones function Jr(q, n), colored by the n-dimensional irreducible representation of SL(2, C); see 
for instance [54]. The HOMFLY polynomial gives rise to an element of ft~ 2 Q[[r, ft 2 ]] using the 
substitution 

t = e T = e Nh , q = e h . (1.2) 

The same substitution works for the colored Jones function of a knot and is the content of the 
so-called loop expansion of the colored Jones function 

oo 

J K (q, n) = J2 FkAtW G Q[[r, h}}. (1.3) 

n=0 

The above loop expansion is really perturbation theory of Chern-Simons theory along an Abelian 
flat connection on the knot complement, see [51]. For a detailed discussion on the meaning of 
the above expansion, see [511 128[ 122] , I21j. Moreover, its planar limit Fk,o(t) can be identified 
with the inverse Alexander polynomial Ak 

Fk,v{t) = -^— y (1.4) 

This is the content of the Melvin-Morton-Rozansky Conjecture, shown in [3]. In addition, for 
every n S N, we have 



where P K ,n{t) G %[t ±l } [Ep. Since A^-(l) = 1, it follows that the radius of convergence of Fk,u{t) 
at t = is positive and independent of n. 



1.3 Power series uniformly Gevrey-1 

In order to state our results, we need to formalize the analytic properties of the free energy of 
a 3-manifold. Recall that a formal power series 

oo 

f(h) = y £a n h n 

n=Q 
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is Gevrey-s if there exists a positive constant C > such that 



On I < C n n\ s 



for all n > 0. 



Definition 1.1. Consider a formal power series 



oo 



f(x,e)= 5„(x)e n Ge- 1 C[[x,e]] 



(1.5) 



n=-l 



of two variables (x,e). We say that f(x,e) is Gevrey-1 with respect to e, uniformly with respect 
to x = (in short, (x, e) Gevrey-1) if there exists a constant C > so that 



Examples of power series f(x,e) (x,e) Gevrey-1 are the WKB solutions of difference or 
differential equations with a small parameter; see for example [1], where the authors call such 
series uniformly pre-Borel summable. The loop expansion of a knot is (r, K) Gevrey-1; see |22j . 

Observe that if a power series f(x, e) given in (jl,5p is (cc, e) Gevrey-1, then for all n € N, the 
formal power series S n (x) is analytic in a common neighborhood of x = 0, independent of n. 

1.4 Statement of our results 

For a closed 3-manifold M, let Zm denote the LMO invariant of M, [H], which is an infinite 
power series of vertex-oriented trivalent graphs, and represents perturbation theory along the 
trivial flat connection of a 3-manifold. Given a metrized Lie algebra q, one can replace graphs 
by rational numbers, keeping track of their number of vertices, and thus create a power series 
(W s o Zm)(K) £ Q[[^]]- Let Zm(N, h) denote the LMO invariant of M, composed with the $l N 
weight system; see [3T]. It is easy to see that Zm(N, K) is a formal power series in N and h with 
constant term 1, and its logarithm 



which is by definition the free energy of M, can be written in the form (jl.ip . where r = iVft. 

Theorem 1.1. For every closed 3-manifold M, the free energy Fm(t, K) is (r, h 2 ) Gevrey-1. 

Theorem 11.11 follows from combining a presentation for the LMO invariant given in Theo- 
rem [L2] with the definition of the gtjy weight system, together with a crucial estimate on the 
number of rooted maps in arbitrary genus obtained by Bender-Gao-Richmond following work 
of Goulden-Jackson; see [TJ [32] and Corollary 14.11 below. 

Theorem 1.2. For every closed 3-manifold M , we can write its LMO invariant Zm in the form 



\[x k ]S n {x)\ <C n+k n\ 



(1.6) 



for all n, k G N. Here, [x k ]g(x) denotes the coefficient of x k 



in a power series g(x). 



F m (t,H) = log Z M (N,H) 



(1.7) 




(1.8) 



r 



where we are summing over the set of trivalent graphs T and 



or I < C 



'■n 



(1.9) 



M 



for all T of degree n, where Cm > is a constant that depends on M . 
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Theorem 11.11 has the following corollary proving 't Hooft's conjecture. 

Corollary 1.1. For every closed 3-manifold M , the power series Fm, 9 (t) are analytic in a com- 
mon neighborhood of t = 0, independent of g. 

We now explain what happens when we specialize N to be a fixed natural number, e.g. N = 2. 
In [20], the first two authors proved the following theorem for a fixed metrized Lie algebra q. 
Let Zm denote the LMO invariant of M and let W s denote the corresponding weight system. 

Then, (W a oZ M )(h)eq[[h}}. 

Theorem 1.3 (|20|. Theorem 3). For every metrized Lie algebra and every rational homology 
sphere M, (W^o Zm){K) is Gevrey-1. 

Corollary 1.2. Theorem W.W implies Theorem 11.31 for q = q[ n for every fixed N. 



1.5 Some calculations 

As a concrete illustration of our results, we can compute the free energy of a Lens space. Let 
L(d, b) denote the Lens space obtained by d/b 6 Q surgery on the unknot in S 3 . 
Recall the a-polylogarithm function 

oo n 
n=l 



for a € R, defined by the absolutely convergent series for \x\ < 1 and analytically continued in 
C \ {0, 1}. For a detailed discussion, see [48] and also [14]. Let B n denote the nth Bernoulli 
number defined by the generating series 



oo R 
Or, 



1 ^ n! 



n=0 



Theorem 1.4. Consider a Lens space M = L(d,b). Its free energy is given by 

F M , g (r) = (2 ff -l)^ [ (d 2 "^Li3^(e^)-Li3_ 25 (eH)+a,(r), (1.11) 

where 



r 2 . . .., r 3 



2 log d - (d 2 - 1)C(3) + \ Lm — if g = 0, 

^(l-0 + ^logci-A £(dj6) I if g = l, ( L12 ) 
if g > 2. 



It follows that for all g > 2 (resp. g = 0, 1), i*M,g has analytic continuation as a meromorphic 
(resp. multivalued analytic) function on C \ d %l^ (resp. C \ ^hr ), where 



Z(l) = 27r»Z, Z*(l) =Z(1)\{0}. 



(1.13) 
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2 The LMO invariant and its i°° Gromov norm 
2.1 The £°° Gromov norm 

The purpose of this section is to prove Theorem 11.21 extending the Gromov norm techniques 
from [2D]. We will assume some familiarity with [20J. To begin with, the LMO invariant Zm 
of a closed 3-manifold takes value in a completed graded vector space -4(0) of vertex-oriented 
trivalent graphs, modulo some linear homogeneous AS and IHX relations, where the degree of 
a graph is half the number of vertices. The LMO invariant gives a meaningful definition to 
the Chern-Simons perturbation theory along a trivial flat connection, and the trivalent graphs 
mentioned above as the Feynmann diagrams of a 3 gauge theory with Chern-Simons action. 
For a detailed discussion, see [39] . The linear AS and IHX relations are the diagrammatic version 
of the antisymmetry and the Jacobi identity of the Lie algebra of the gauge theory. 

Let -4 n (0) denote the subspace of -4(0) of degree re. Then, -4 n (0) is a finite dimensional 
vector space spanned by the finite set of vertex-oriented trivalent graphs with 2re vertices. Of 
course, this spanning set is not a basis, due to the linear AS and IHX relations. This motivates 
the following concept of the £ p Gromov norm, extending the one in [201 Definition 1.1]. 

Definition 2.1. Consider a vector space V and a spanning set b. Fix p G [1, oo]. For v G V, 
define the £ p norm by: 



mf (£,-N P ) 1/P if pe[l,oo) 
inf maxj \cj\ if p = oo, 



(2.1) 



where the infimum is taken over all presentations of the form v = £ ■ CjVj, Vj G b. 

We now apply the above definition to each of -4 n (0) and combine them into a single power 
series. 

Definition 2.2. 

(a) Fix p G [l,oo] and re G N and consider the vector space -4 n (0) spanned by the set b of 
oriented trivalent graphs of degree n. For v G -4 n (0), we denote by \v\ p the norm of v. 

(b) If v G -4(0), we define 

oo 

\v\ P (h) = ^\n n (v)\ p h n €Q[[h]] (2.2) 
n=0 

where ir n : -4(0) — ► -4 n (0) denotes the projection in -4 n (0). 

(c) We say that v G -4(0) has Gevrey-s p-norm if |f|p(fr) is Gevrey-s. 

2.2 A brief review of [20] 

The main result of [20] was the following theorem. 



Theorem 2.1 (|20j. Theorem 2). For every integral homology sphere M, the power series 
\Z M \i{h) G Q[[h}] is Gevrey-l. 

Theorem l2.1l is a convergence property of a universal finite type invariant of homology spheres, 
and uses a well-known universal finite type invariant of links in S 3 , namely the Kontsevich 
integral. For an introduction to the notion of finite type invariants of knots and 3-manifolds, 
see [2l|50]. Since our proof of Theorem 11.21 will use some of the ideas of the proof of Theorem 12. 1[ 
we begin by recalling the proof of Theorem I2.1[ which goes as follows. First, we begin with the 
following description of the LMO invariant from [5] . 
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• Fix a presentation of an integral homology sphere M as surgery on a unit-framed boundary 
link L in S s . 

• Choose a presentation of L as the closure of a framed string link T. 

• Consider the normalized Kontsevich integral Zt, which takes values in the completed Q- 
vector space A(*x) of vertex-oriented Jacobi diagrams with legs colored by a set X in 1-1 
correspondence with the components of T. 

• Perform formal diagrammatic Gaussian integration 



to Zt- I.e., separate out the strut part of Zt (which is an exponential of Jacobi diagrams 
with no trivalent vertices) and join the legs of the remaining diagrams using the inverse 
linking matrix of the strut part. 

• Finally, define 



where denotes the ±1 framed unknot in S 3 and o~±(L) denotes the number of positive 
(resp. negative) eigenvalues of the linking matrix of L. 

Next, we extend our notion of £ p norm to all intermediate spaces of Jacobi diagrams (with 
or without skeleton, and with or without symmetrized legs). 
Next, we show that the Kontsevich integral is Gevrey-0. 

Theorem 2.2 ([20], Theorem 8). For every framed tangle T, \Zt\i(K) G R[[^]] is Gevrey-0. 

This follows from the definition of the Kontsevich integral using the KZ associator, together 
with the following lemma. 



Lemma 2.1 ([20], Proposition 2.6). IfU denotes the zero-framed unknot, then \Zu\x(fi) £ K[[^]] 

is Gevrey-0. 

Finally we use the following key lemma. 

Lemma 2.2 ([20], Lemma 2.10). For every boundary string link T with framing ±1, 
|J Z T dX\i(K) E Rp]] is Gevrey-l. 

The condition on the string link being a boundary one is to ensure that the strutless part 
of Zt contains no trees; see [33]. Lemma [2.2l takes care of the normalization of the LMO invariant 
given by (]2.3p . Ignoring technicalities, the main point in the proof of Lemma 12. 2 1 is the following. 

If G is uni-trivalent graph of degree n + k (i.e., with 2n + 2k vertices) with 2k legs, then 

there are (2k — 1)!! = 1-3-5 (2k — 1) ways to pair the legs of G together. If G has no tree 

components, then k <n thus 

(2k - 1)!! < (2n- 1)!! < C n n\. 

Remark 2.1. If we replace | • |i by | • | p for p £ [l,oo), Theorems 12.11 [2~2l and Lemmas 12.11 12.21 
remain true. 

In the next section we discuss what happens when we use the | • |oo rather than the | • |i norm. 




J Z L dX 



(2.3) 
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2.3 Proof of Theorem Q 

We now show how to modify the statements of the previous section in order to show the following 
reformulation of Theorem 11.21 In that theorem we need to fix a closed 3-manifold. When M 
has positive Betti number, the LMO invariant simplifies a lot, and it is possible to give a direct 
proof by an explicit formula. For example, when the first Betti number &i(M) is greater than 3, 
then the LMO invariant vanishes and Theorem 11.21 obvious holds. When 1 < b\(M) < 3, the 
LMO invariant can be directly computed, see [34} 150]. When b±(M) = 1, the LMO invariant 
is computed by [24] in terms of the Alexander polynomial of M , and Theorem 11.21 follows 
immediately. Thus, it suffices to assume that M is a rational homology sphere, and further 
(using the multiplicative behavior of the LMO with respect to connected sums), it suffices to 
assume that M is obtained by surgery on a rationally framed boundary link in S 5 . In that case, 
the LMO invariant can be computed via the Aarhus integral. 

Theorem 2.3. For every rational homology sphere M, |^Af|oo(^) £ ^[[^]] Gevrey-0. 

This theorem follows from the following. 
Theorem 2.4. For every framed tangle T, \Zt\oo{K) G ^[[^]] is Gevrey-0. 

This follows easily from Theorem 12.21 and the fact that |t>|oo < |t>|i for all v S A n ($). 

Lemma 2.3. For every framed string link T with invertible linking matrix \ j Zr<iX| 00 (^) £ 
R[[h]] is Gevrey-0. 

The proof of this lemma is as follows. First we write Zt as follows, using Theorem 12.41 



summation is over the set of uni-trivalent graphs with no strut components and legs colored by 
the components of T and 



for some constant C > 0. If G is a trivalent graph of degree n, it has 3n edges. If G comes from 
pairing the legs of a uni-trivalent graph T, then T is obtained from G by cutting some (let's say 
k < 3n) of the edges of G in half. Such a uni-trivalent graph T has degree n + k (i.e., 2n + 2k 
vertices) and 2k legs and its coefficient in Zt is bounded by C n+k < C 4n . Since G has 3n edges, 
there are at most 2 3n such graphs T. Thus, pairing the legs of the strutless part of Zt, we obtain 
that 



Zt = exp I - 




where (Zy) is the linking matrix of T, 



is a strut colored by the components i and j of T, the 



cr| < C de ^ 




G 



where 



cg\ < O 



</dcg(G) 



for some C > 0. This concludes the proof of Theorem 11.21 
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3 The gi N weight system 



The LMO invariant of a closed 3-manifold takes values in a completed vector space A($) spanned 
by trivalent graphs. To get a numerical invariant with values in Q[[/i]], we need to replace every 
graph by a combinatorial weight. This is exactly, what a weight system does. More precisely, 
given a Lie algebra g with an invariant inner product, there is a weight system Q-linear map: 



W s : A((l>) — Q[[h}} 
discussed at great length in [2]. In [2], the following description of W s is given. 

W glN :A(®)^®[[N,h}} 
is defined by: 



(3.1) 



M 



where 



the sum is over all possible markings M of the trivalent vertices of D by or 1, 
sm is the sum, over the set of trivalent vertices, of the values of M, 

S_d,m denotes the X-marked surface obtained by thickening the trivalent vertices of D as 
follows: 







and thickening the edges of D, and connecting up into to a surface. It turns out that £d,m 
is well-defined, oriented and perhaps disconnected, 

• 9d,m and ^d,m denote the genus (i.e., sum of the genera of its connected components) and 
the number of boundary components of T,(D,M), 

• it turns out that the degree of D is Ig^^i — 2 + 6d,a/- 
For example, we have 

^ 0[jv (9) = 2(iV 3 -Ar)/ l , 

where is the obvious planar trivalent graph with counterclockwise cyclic order at each vertex. 
The next lemma follows directly from the above definition of the W g i N weight system. 

Lemma 3.1. IfT is a connected trivalent graph of degree n, then 

W slN (T)= Pr (N)h n , (3.2) 

where pr{N) £ Z[N] is a polynomial in N of degree at most n + 2 and l l -norm at most 2 n . 
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4 Asymptotics of the number of rooted maps of arbitrary genus 

In a series of papers in various collaborations, Bender et al studied the exact and asymptotic 
number of rooted maps on a surface; see [HE]- Recall that a map (G, S) is a graph G embedded 
in a connected, oriented, closed topological surface S in such a way that each maximal connected 
component of S \ G is a topological disk. A map is rooted if an edge, a direction along the edge, 
and a side of the edge are distinguished. Let T g {n) denote the number of n-edged rooted maps on 
a surface of genus g. In [6], Bender-Canfield, using generating functions, and the Darboux-Polya 
enumeration method, give the following asymptotic expansion of T g (n). 

Theorem 4.1 ([6]). For g fixed and n — > oo we have 

T g (n) ^tgn 5 ^ 1 )/ 2 ^, (4.1) 

where the constants t g are computable via computable non-linear recursions. In particular, 
2 1 7 

^0 — —}=, t\ — — , t2 



/F 24' 4320^' 

Many interesting families of maps satisfy asymptotic formulas of the form 

at g {f3nf^)n^_ 

This is discussed extensively for example in [551 Section 4]; see also [29|. I30J . For over 20 years, the 
constants t g remained hard to compute or estimate, partly due to the complexity of their non- 
linear recursion. A breakthrough was achieved recently. Using work of Goulden-Jackson [32] 
on the KP hierarchy, Bender-Gao-Richmond simplified the nondinear recursion relation for t g , 
and obtained the following. 

Theorem 4.2 ([7]). We have 

40sin(7r/5)if / 1440<? \ ~ 9/2 

ta {—) ' ( } 

where K = 0.10486898772254091800 ... is a constant. 

The numerical value of K was given incorrectly at [7j- An exact value of K is given in 
equation (|A.16p in Section IA1 below. For our purposes, it suffices to note that a trivalent graph 
of degree n has 3n edges. Taking into account the choice of a root, it follows that 

Corollary 4.1. There are at most 12nT g (3n) connected trivalent graphs of genus g, where T g (n) 
satisfies (j4TT|) and P~2|) . 



Remark 4.1. Theorem 14.11 has an analogous statement (with a different constant p g instead 
of t g ) for rooted graphs in unoriented surfaces; see [6j Theorem 1]. Without doubt, there is an 
asymptotic expansion for p g analogous to t g . In fact, one can write down explicitly a non- linear 
recursion relation that p g ought to satisfy, and study its exact asymptotic behavior. This will 
be discussed in a separate publication; see [27] . 



5 Proofs 

5.1 Proof of Theorem ll.il 

Let Zm(N, K) denote the LMO invariant of M, composed with the weight system of gl^r. Then, 
we can write 

log Z M (N,h)= a M , g ,dN d h 2 ^ 2+d . (5.1) 

2g-2+d>0,d>0 
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Recall that the free energy of M is defined by (jl.7p 

f m (t, h) = log z m (n, h), 

where r = Nh. It follows that the free energy has the form (jl.ip where 

Fm, b {t)= Yl «%/eQ[Hj. (5.2) 

d:2g-2+d>0,d>0 

On the other hand, Theorem 11.21 implies that we can write 

log z M = c r • r, 

r 

where the sum is over a set of connected trivalent graphs, where there exists a constant C such 
that 

|cr| < C n (5.3) 
if r has degree re. Applying the W g i weight system, and using Lemma 13. II it follows that 

oo 

F m (t, n) = ^^c r -p r (iV)/i n , 
n=i r 

where the T summation is over a set of connected trivalent graphs of degree n. Thus, 

a M , g4 = Y,cr-lN d } Pr (N), (5.4) 
r 

where the sum is over the set of connected trivalent graphs T of degree re = 2g — 2 + d that 
embed on a surface of genus g. Lemma 13.11 estimates the coefficients [N d ]pr(N) and since 
2g — 2 + d = re, it implies that 



\aM, g ,d\ < C™ ^ l c r| 



r 

for some constant C\ > 0. Equation (|5.3p together with 2g — 2 + d = n implies that 

\a>M,g,d\ < C% ^2 1 
r 

for some constant C2 > 0. The above sum is over the set of trivalent graphs of degree n = 
2g — 2 + d (and thus, 3n edges) that embed on a surface of genus g. Corollary 14. II implies that 

(3n) 5 \ 9/ * (2g-2 + df^ 



J2 1 < 8nT g (3n) < C{ 1 ( - 

T V 



< r n — 

9 J - 2 



(2g-2 + df^ (2g)^ 
2 (2g) 5 9/ 2 g9/ 2 ' 

Since (1 + a/n) n < e a , we have 

{2g-2 + dfy' 2 5(d _ 2)/4 
{2gf9/2 

Moreover, 

Since re = 2g — 2 + d, it follows that 

\a M , g ,d\ <Cl +d (2g)l 
and concludes the proof of Theorem 11.11 
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5.2 Proof of Corollary ITT21 

Fix a natural number iV£N. It suffices to prove that (W g i N o log Zm)(K) G Q[[^]] is Gevrey-1. 
Equation (15. 1ft implies that 



m(W glN o logZ M )(ft) = J2 a M , 9>d N d . (5.5) 

g,d:2g— 2+d=n 

The proof of Theorem 11.11 implies that there exists C > so that 

\a M , g ,d\ <C 9+d (2g)l. 
Since 2 5 - 2 + d = n, C^ d < C* n/2 = C% and (2g)\ < C$n\. Thus, 

|coeff(/f\ {W g[N o logZ M )(^))| < C>! iV d < C 4 n n!AT n+2 (n + 2). 

g,d:2g— 2+d=n 

Since N is fixed, it follows that 

\[h n ](W slN o log Z M )(h)\ <C n n\ 
which concludes the proof of Corollary 11.21 

6 The free energy of a Lens space 

In this section we will prove Theorem 11.41 The next proposition computes the image of the 
LMO invariant of a Lens space L(d,b) under the weight system W s . 

Let g denote a metrized Lie algebra with inner product (•, •) and let denote the positive 
roots of g, | $-|-| denote the cardinality of <&+ and p denote half the sum of the positive roots 
of 0. Let c denote the product of the quadratic Casimir of with the dimension of 0. Let Am 
denote Casson invariant. 

Proposition 6.1. With the above assumptions we have 

iWa ZLmm . exp (^v) j + «f e , 6 ,) 

Proof. There are two proofs of this proposition. One proof uses 

(a) The existence of the Ohtsuki series (which come from the Reshetikhin-Tuarev invariants 
of 3-manifolds). This was shown by Ohtsuki in [49] for PSU(2), and by the second author 
in [12] for PSU(iV) and more generally in [53] for all simple Lie algebras. 

(b) The computation of the Ohtsuki series for Lens spaces, done by Takata in [52j . 

(c) The identification of the left hand side of (|6.ip with the Ohtsuki series. This was the 
result of the unpublished fourth part of [5] and independently an unpublished work of the 
second author that was recently completed by Kuriya [10], see also |50| . 



An alternative proof uses the computation of the LMO invariant for Lens spaces by Bar-Natan 
and Lawrence [1] . Let us give the details of this proof. We will use the notation from [1] . In [H 
Proposition 5.1] it is shown that the LMO invariant of L(d,b) is given by: 

S(b/d) \ (Qx,Q x /d)x 



Z L (d,b) = exp ^ 4g 8 j 



(ri x , Qx) x 
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Moreover, —S(b/d) = 12\ L (p d y Now, we apply the weight system W g . We have 

w s (&)(h) = c h, 

n ^ 

a£$ + \ if/ i 

and 

w B ((n x ,n x/d ) x )(h) = w s ((n x ,n x ) x )(h/d). 

The result follows. ■ 
6.1 Some special functions 

Let us introduce some auxiliary functions which have already appeared in the LMO invariant of 
the Lens spaces and which are important ingredients of the proof of Theorem II, 41 Consider the 
functions / and / defined by 

~ sinh(x/2) ~ 

f( x ) : = 77^ — . and f(x) := log/ (x). 

x/2 

It is known that 

00 2 4 fi 



^ 24 2880 181440 

k=l 

where &2fc's are the modified Bernoulli numbers and are related to the ordinary Bernoulli num- 
bers -Bfc as follows: 

2k ~ 2k(2k)V 
The first polylogarithm is given by 

Lix(x) = -log(l -x). (6.2) 
The relation between f(x) and Lii is the following 

/(x) = -Li 1 (e a; )-|-log(-x). (6.3) 

The left hand side is an analytic function in the unit disk |x| < 1 with an analytic continuation 
as a multivalued analytic function in C \ 2ni(Z \ {0}). The right hand side is well-defined in the 
unit disk minus [0, 1), and extends as a continuous function over the cut [0, 1). Then, both sides 
agree on the unit disk \x\ < 1. Equation (|6.3[) follows from the following easy computation 

/V/2 _ e -^/2\ ( g-z/2 \ 
f{x) = log = log (1 - e*) 



and equation (|6.2|) . 

From the power series definition of the polylogarithm, it is easy to see that for all a £ R we 
have 

^U a (e x )=U a _ 1 (e x ). (6.4) 
Moreover, when a is a negative integer, then 

where P a {x) G 7L\x\ is a palindromic polynomial of degree —a. It follows that for a a negative 
integer, the function Li a (e r ) is an even meromorphic function on C with poles at Z(l), where 
the latter set is defined in (|1 . 13|) . 
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6.2 Proof of Theorem Q 

This section is devoted to the proof of Theorem 11.41 Let us observe that the gt^y and the 
sIn weight system agree on all nonempty trivalent graphs (since the structure constants of an 
Abelian Lie algebra vanish, and are placed at the trivalent vertices of a trivalent graph), and 
that the logarithm of the LMO invariant of a closed 3-manifold is a series of nonempty connected 
trivalent graphs. Therefore, to compute the free energy, we can work with the sIat weight system. 
Recall that the set $>+ of positive roots of sIn has N(N — l)/2 elements. From Proposition 16.11 
one has 

4f + /(Kris) 

Taking the logarithm, we get 

F M (T,h) = h^±c BlN + E (f«a,p)h/d)-f((a,p)h)). 

Note that 



When a runs through the set $+, (ct,p) takes the value j the total of N — j times, for every 
j = 1, . . . , N — 1. Hence 



A 7V ~ 1 

F M (r, h) = h + £ ((AT - - (JV - • (6.6) 



Now we have 

JV-l N-l oo oo AT-l 

J2(n- j)f(jh) = j2(N-j)j2 b2kj 2k n 2k = e ^ J2(n- j) j 

3=1 j=l k=l k=l j=l 



2k 



Using the well-known identity expressing the sum of powers by a Bernoulli polynomial 
N-l , k 



o— i o— n V / 



k+l-s 



j=l s=0 

it follows that 

N-l k 



Y (N-i) j 2k = V (2fc)!(1 ~ 2g) B 2o N 2k+2 - 2 9. 
pi f^(2 9 )l(2k + 2-2g)l 29 



Therefore, using N = r/h, we get 

N-l oo 



^ (1 - 2g) B^ 2 -^ (2AQ! 2fc+2 _ 2 , 

^ (2o)! ^ (2fc + 2-2 5 )! 2fc 

9=0 V y; fc=max{ 9 ,l} V y/ 
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3=0 



(l-2g)i? 2 ^ 2 - 2 ^(2; + 2 g )! ^ 2l+2 



(6.7) 



with a minor variation when g = 0. Let us define the auxiliary functions F g by 
' (2Z + 2 5 )! ; 



(2/ + 2)! 



J2g+2V' 



r 2l+2 



(6.8) 



Then, equations ()6.6p . (|6.7f) and (|6.8f) and the replacement of (ft, r) by (h/d,T/d) imply that 



B 



F Mtg (r) = (1 - 2 5 )^ (d 2 ^F g (r/d) - F g (r)) + '^(rH gfi - r5 gA ). 

We claim that the functions F g are given by 

t 2 r 3 St 2 7r 2 r 

- T lo g (-r)-- + --- + C(3) if , = 0, 



F ff (r) = -Li 3 _ 29 (e T ) + < 



- - -log(-r) 
(2 5 -3)!r 2 " 29 



2g~2 



2g-2 



if 9 = 1, 
if q > 2. 



Theorem 11.41 follows easily from equations (|6.9p and f|6. lOj) . 

It remains to prove (|6.10p . Equation (|6.8p and ()6.3p implies that for g = 1 we have 

*i(r) = /(r) = -Lix^) - I - log(-r). 



From equation (|6.8|) . it is easy to see that for 5 > 1 we have 

F g {r)=d^- 2 F 1 {r)-{2g-2)\b 2g . 2 . 
Since for all a € R we have 

<9 T Li Q (e T ) = Li a _i(e r ). 
Equations (|6. 1 1 j) and (|6.12p imply that for g > 2 we have: 

F g {r) = -Li 3 _ 29 (en + (2g - 3)!r 2 ^ - ^*=|, 

where we have used the fact that (2g — 2)!62 9 = B2 g /(2g — 2). 
When 5 = 0, equation (|6.8p implies that 



(6.9) 



(6.10) 



(6.11) 



(6.12) 



(6.13) 



F (t) = d~ 2 f(r) 

up to a linear function of r, where d~ l denotes integration with respect to r. Integrating (|6.13|) 
twice using (16. 4h . and matching the first three coefficients of the Taylor series of both sides at 
r = implies (|6.10p for g = 0. This concludes the proof of Theorem 11.41 

An alternative proof of equation (|6.10|) follows from Gopakumar-Vafa |3H Appendix]. See 
also [13116]. 

Yet another proof of equation (|6.10p follows from |13t Proposition 1]. ■ 



Remark 6.1. The reader may compare Theorem 11.41 with the BPS formulation of [31 j and 
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7 Future directions 

7.1 Analyticity of the free energy of a matrix model 

The notion of free energy exists for U(N) gauge theories where N is arbitrary. Toy models of 
those theories are the so-called matrix models studied by the French school \12\ [TO] . In a future 
publication [2U], we will study the analytic properties of the free energy of a matrix/multi- 
matrix model with arbitrary potential. On the other hand, it was shown in [25] that given 
a closed 3-manifold M presented by surgery on a framed link L in S 3 , there exists a multi-trace 
potential Vl whose free energy coincides with the free energy of M. This observation gives 
geometric examples of different potentials with equal free energies. Combined with the future 
work of [26] , this may give another proof of our main Theorem 11.11 

7.2 Analytic continuation of the free energy in the complex plane 

This section is motivated by the questions that were raised in [23]. Fix a closed 3-manifold M, 
and let pM,g denote the radius of convergence of the power series Fvi !ff (r) at r = 0. Theorem ll.il 
implies that ixii{pM,g | g > 0} > 0. 

Question 7.1. Fix a closed 3-manifold M. Is it true that the radius of convergence of Fm^{t) 
at t = is independent of g? Do the power series Fm^t) admit analytic continuation as 
multivalued analytic functions in the complex plane minus a set of singularities, independent 
of 9? 

Notice that if we replace the closed 3-manifold M by a knot K, and consider its free energy 
from the loop expansion of K, then the answer to this question is positive. This was discussed 
in Section fl. 21 In that case, Fr^t) are rational functions of e T and have analytic continuation 
in the complex plane minus the logarithms of the roots of the Alexander polynomial of K. 

Question 7.2. Can you compute the free energy (or even its planar limit) for any closed hy- 
perbolic manifold? 

7.3 The double scaling limit 

In the physics literature, it is customary to expect that F g (r) has an expansion around a g- 
independent singularity tq as follows: 

F g {r) = (r - T yv +S c g (2g)\ + 0((r - r )™ +1 ). 

In that case, one considers the single-variable power series 

oo 

f{x) = ^T t c g (2g)\x9, 

9=0 

which is called the double scaling limit of F. In the case of matrix models, 2-dimensional gravity 
and random matrices, it turns out that f(x) satisfies non-linear differential equations which 
are a specialization of the KP hierarchy. For a lengthy discussion, see [16], [551 Section 4] and 
references therein. 

Question 7.3. Does the double scaling limit of the free energy of a closed 3-manifold satisfy 
a specialization of the KP hierarchy ? 
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7.4 The O(N) and Sp(N) theories 

In the present paper, we defined the free energy of a closed 3-manifold using the U(N) gauge 
theory. We could have used the O(N) or the Sp(N) gauge theory. In that case, the weight 
systems W 0N and W s p N lead to trivalent graphs that embed to unoriented surfaces, see [2]. 
Theorem 11.11 holds in that case, assuming an asymptotic formula for the number of rooted 
unoriented maps; see Remark 14.11 A precise asymptotic expansion for the constants p g of 
Remark |4. II will be discussed in a separate publication [27] , 



A Asymptotics of t g and Painleve I 

In this Appendix we show that the recursion relation for t g found in [7J is closely related to 
an asymptotic, formal solution to the Painleve I differential equation. This makes possible to 
refine the asymptotic behavior obtained in [7J. 

We first recall some results from [7J . The asymptotics of t g is obtained from the asymptotics 
of two auxiliary sequences. The first one is f g , which is defined by the recursion 



/fi 

fg = - 4) (5 5 - 6)/ g _! + eVE^fhfg-h, (A.l) 

h=X 



and the initial condition 

A - (A.2) 
The sequence t g is related to f g by 

f g = 24" 3 /V/ 2 r (^J 1 ) t B . (A.3) 
Moreover, [7J introduces another sequence u g , defined by 

Using (|A.ip . it is shown that u g approaches a constant K as g — > oo, and this leads to the 
main result of [7J concerning the asymptotics of t g . The constant K is in principle only known 
numerically. 

To start our analysis, we define yet another sequence a g by the relation 



V6 
2 



f g = ibC%, C = V> (A.5) 



so that the recursion (|A.1|) becomes 

25(g-l) 2 -l 1^ 1 ,. n 

ag = 48 a H _ 2L fl ' fl H' ao = l. (A.6) 

We now consider the formal power series 

oo 

Mz)=Z^ a 3 Z ~ 59/2 - ( AJ ) 
9=0 
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It is then easy to see that the recursion (|A.6P implies that (fio(z) satisfies the following differential 
equation 

f ~ If" = z, (A.8) 

which is the well-known Painleve I equation. Therefore, the differential equation presented at 
the end of [7] is Painleve I in disguise. 

We can now use properties of Painleve I to give a precise description of the asymptotics 
of a g , and in turn of f g and t g . This follows from the general results of [15] together with the 
computations of [18] for Painleve I; see also |38^ Section 3]. We will then content ourselves with 
a statement of the asymptotic behavior of a g and a sketch of the main ideas. 

The coefficient a„ has the following asymptotic behavior as g — ► oo: 




r 2g-± S 1 + V , - ■ (A. 




7T 

In this expression, 

^4 = ^, * = (A.10) 
5 lyj n 

and the \i\ are defined by the recursion relation 

& = T7T7f7 \ 2J WV-k+W -( l -Yo)( l - Jo) J , Mo = 1- (ATI) 

It is understood here that a n / 2 = if n is not an even integer, otherwise the coefficient a n / 2 is 
given by the recursion (|A.6p . The quantities S and fi[ have a very clear interpretation in terms 
of the differential equation (|A.8[) . It is known (see for example [18] ) that the Painleve I equation 
exhibits a non- linear version of the Stokes phenomenon across the line Argz = 0. S turns out 
to be the corresponding Stokes parameter. The exact value of S was first obtained in [37] using 
the Riemann-Hilbert approach, further justified by [19]. A calculation of the Stokes parameter 
using a Borel summation approach and WKB was given by [53] . The coefficients \i\ arise as 
follows. The Painleve I equation admits a so-called formal trans-series solution of the form 

00 

4>(z) = fa{z) + z$J2?M*), £ = Cz~ 5 / 8 e- Az5/ \ (A.12) 

where <Pq{z) is the formal power series ()A.7p . C is a constant and (j>g(z) are formal power series 
of the form 

oo 

= X>,k z " 5 * /4 - (A.13) 

fc=0 

normalized by (j>i(0) = 1. The Painleve I equation gives a series of recursion relations for the 
coefficients (f>g & in terms of a n and 4>£> t k' with £' < I. The coefficients of the first series in the 
trans-series, (f>i,ki satisfy precisely the recursion (jA.llh . therefore fii = (f>\ t i. One finds, for the 
very first terms, 

, f \ 5 _ 5 75 _s 341329 _is 

6i(z) = 1 =z * -| z 2 _2i 4 -| . (A. 14) 

U 7 64^ 8192 23592960^3 
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The fact that the value of the Stokes parameter S together with the coefficients of the first 
trans-series correction (j>\{z) lead to an asymptotic behavior of the form ()A.9|) has been known 
for a long time in the physics literature on large-order behavior in perturbation theory; see for 
example [44;. A rigorous proof for the case of Painleve I was given in [38[ Section 3]. For a 
general approach to this problem for a class of generic non-linear ODEs, see [15]. Using these 
results we can now revisit the asymptotics of the sequences u g and t g . One finds indeed that 



as g — ► oo. The constant K can now be determined exactly from the asymptotics of a g . Its 
value is, 



In fact, the numerical approximation to the value of K presented in [7J is not very precise, and 
can be substantially improved by using standard techniques like Richardson transforms. Of 
course, the results above give all the subleading corrections to the asymptotics of u g (and t g ) in 
powers of 1/g. 

Let us use the Painleve equation to answer a question posed by [7J concerning the sequence t g . 
The question asks whether the sequence (t g ) is holonomic, i.e., satisfies a linear recursion in g 
with coefficients polynomials in g. It is easy to show that if t g is holonomic, then the generating 
series 4>o(z) satisfies a linear ODE with coefficients polynomials in z. It follows that (f>o(z) (or 
more precisely, its Borel transform in a suitable ray) has analytic continuation as a multivalued 
analytic function in C minus a finite set of singularities. On the other hand, it is known that 
any (truncated) solution of Painleve I is meromorphic (see for example [111 [36] ) , with infinitely 
many singularities that are asymptotic to one of the rays arg(z) = 2irk/5 for k = 1, . . . , 5 (see 
for example, [PL] [19]). 

Corollary A.l. The sequence (t g ) is not holonomic. 

We end this Appendix by pointing out some connections between the results of [7J [32] and 
the physics literature on matrix models and two-dimensional gravity. It has been known since 
the pioneering work of p21 [10] that matrix models are a useful tool to count maps of genus g and 
with n vertices. The resulting generating functions at fixed genus g are all analytic functions in 
the vertex-counting parameter, and they have a finite radius of convergence which is common to 
all of them. The sequence f g in [7J is given indeed by the coefficients of the leading singularities 
for one of these generating functions. It is known however from the physics literature on the 
so-called "double scaling limit" that these coefficients are governed by the Painleve I equation 
(see j!6j). and in fact, as we have seen, this is the differential equation governing the sequence f g . 
A study of the asymptotics of various generating functions appearing in matrix models, including 
the Painleve I equation, can be found in [47] . 
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